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ABSTRACT 

We report the discovery of the near-infrared and optical afterglow of the short-duration gamma-ray burst 
GRB 070724A. The afterglow is detected in UHK^ observations starting 2.3 hr after the burst with Kg = 19.59 ± 
0.16 mag and / = 23.79 ±0.07 mag, but is absent in images obtained 1.3 years later Fading is also detected in 
the /Ts-band between 2.8 and 3.7 hr at a 4a significance level. The optic al/near-IR spectral index, /3o,nir ~ -2, 
is much redder than expected in the standard afterglow model, pointing to either significant dust extinction, 
Ay°^' « 2 mag, or a non-afterglow origin for the near-IR emission. The case for extinction is supported by a 
shallow optical to X-ray spectral index, consistent with the definition for "dark bursts", and a normal near- 
IR to X-ray spectral index. Moreover, a comparison to the optical discovery magnitudes of all short GRBs 
with optical afterglows indicates that the near-IR counterpart of GRB 070724A is one of the brightest to date, 
while its observed optical emission is one of the faintest. In the context of a non-afterglow origin, the near-IR 
emission may be dominated by a mini-supernova, leading to an estimated ejected mass of 10-* Mq and 
a radioactive energy release efficiency of / ^ 5 x lO""* (for v ^ 0.3c). However, the mini-SN model predicts a 
spectral peak in the UV rather than near-IR, suggesting that this is either not the correct interpretation or that 
the mini-SN models need to be revised. Finally, the afterglow coincides with a star forming galaxy at z = 0.457, 
previously identified as the host based on its coincidence with the X-ray afterglow position 2" radius). Our 
discovery of the optical/near-IR afterglow makes this association secure, and furthermore localizes the burst to 
the outskirts of the galaxy, with an offset of 4.8 ±0.1 kpc relative to the host center. At such a large offset, the 
possible large extinction points to a dusty environment local to the burst and rules out a halo or intergalactic 
origin. 

Subject headings: gamma-rays ;bursts 
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1. introduction 

The determination of sub-arcsecond positions for short- 
duration gamma-ray bursts (GRBs) is of the utmost 
importance for our growing understanding of their red- 
shift distribution, energy scale, host galaxies, and local 
environments. Such localizations require the detection of 
optical, near-infrared, or radio afterglows; or alternatively an 
X-ray detection with the Chandra X-ray Observatory. As 
of August 2009, only 15 short GRBs have been precisely 
localize d in this man n er, an d all were d e tected i n the optical 
band^ (iBerger et al I l200l lFoxetal.1 l200l iHiorth etal.1 



2005bt iLe vanetal.' '2006'; 'de Ugarte Postigo et al.' '200^ 



Roming et a l. 2006; Soderberg et al. 2006; Perlev et al. 2007] 
Stratta et al . 2007; Malesani et al. 2008; Piranomonte et all 

2009: 



2008r IBerger & Kelson 



.Cenko et al., ,2009a 



D'Avanzo et al. 2009; Gr aham et al.ll2009HKuin & HoverstenI 
2009i:.Perley et al., .2009b). 



Follow-up observations of some of these bursts have led to 
the detection and characterization of host galaxies, most of 
them star for ming, and a sm all proportion 20%) non-star 
forming (e.g.. lBergerll2009l) . The precise positions also allow 
us to study the local environment of the bursts within their 
hosts, and current observations point to offsets of ^ 1 ~ 15 
kpc (e .g., Berger et al. 2005; F ox et"ani2005l;ISoderberg et all 
I2006t iD'Avanzo et al.,,2009 ). However, since these are pro- 
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jected positions, and since we generally lack detailed af- 
terglow observations that can shed light on the circumburst 
environment, little is known about whether the bursts orig- 
inate within the inner or halo regions of their hosts (e.g., 
Soderberg et al. 2006; Levesque et al. 2009). 

Rapid optical observations have also been used to place lim- 
its on emission from radioactive material synthesized in a pu- 
tative sub-relativistic outflow associated with a compact ob- 
j ect binary m erger, a so-called Li-Paczynski mini-supernova 
dLi & Paczv riski 1998). Such emission is theorized to have a 
typical peak time of ~ 1 d, and a peak luminosity of ~ 10"*^ 
erg s"', corresponding to m^B ^21 mag at z ^ 0.5. No 
such emission has been d etected to date (e.g., iHiorth et al.l 
l2005al; iBloom et"aLll2006l) . Similarly, no late-time emission 
from Type Ib/c supernova associations have been detected 
(e^ lHjorth etani2005bl; iBloom et al.lf2006l; ISoderberg etaP 
l2006b. 

Most recently, near-IR and optical non-detections of 
GRB 070724A have been used to place limits on emis- 
sion from a putativ e mini-SN associated with this burst 
dKocevski et al.l l2009). Here we report the detection of near- 
IR and optical counterparts of GRB 070724A about 2.3 hr af- 
ter the burst, and show that the afterglow is actually one of 
the brightest near-IR short GRB afterglows detected to date, 
but is one of the faintest in the optical. We use the observed 
fluxes and the unusually red color to investigate the properties 
of the afterglow and/or mini-SN, and to precisely measure the 
location of the burst relative to its host galaxy. Our discovery 
of the afterglow of GRB 070724A suggests that recovery of a 
substantial fraction of short GRB optical/near-IR afterglows 
requires observations to m^B ^ 24 mag within about 0.5 d. 

2. GRB070724A 
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GRB 070724A was discovered by the Swift satellite on 

2 007 July 24 at 10:53:5 UT with a duration of 0.40 ± 0.04 
s jZiaeepour et al.ll2007h . The X-ray afterglow was detected 
with the on-board X-ray Telescope (XRT) beginning 72 s 
after the burst, while no counterpart w as detected with the 
UV/Optical Telescope (IZiaeepour et al.ll2007 ). The X-ray po- 
sition was subsequently determined to a precision of 1.7" ra- 
dius (90% containment). An apparently extended source was 
detected in coincidence with the XRT error circle in Digitized 
Sky Survey images, and tentatively proposed as a possible 
host (Bloom 2007). Subsequent near-IR and optical observa- 
tions from UKIRT, Gemini-North, the Palomar 60-inch tele- 
scope, and the VLT revealed that the source was in deed an ex- 
tended galaxy, but did not uncover an afterglow ( Levan et alj 
120071: ICenko et al.ll2007l: ICovi no et ar'2007). A recent anal- 
ysis using image subtraction on UKIRT, NOT, CTIO 1.3-m, 
Keck, and VLT data similarly reveals no afterglow to limits 
of Fi,,K ^ 30 and /y.i 5, 0-3 AtJy at 3.2 and 22.2 hr after 
the burst, respectively (iKocevski et al.ll2009h . 

Spectroscopy of the galaxy within the XRT error circle re- 
vealed that it is located at a redshift of z = 0.4571, is undergo- 
ing active star formation at a rate of 2.5 M0 yr~', has a lumi- 
nosi ty of Lb » 1 .4 L*, and a metalhcity o f 12 + l og(0/H) « 
8.9 ( lBergerfl2009t see also iKocevski et all ( l2009l) for similar 
results). The resulting isotropic 7-ray energy in the observed 
15-150 keV range is E^j^o ~ 1 .6 x lO'*^ erg. 

3. DISCOVERY OF THE NEAR-IR AND OPTICAL AFTERGLOW 

We observed the field centered on GRB 070724A with the 
Near Infra-Red Imager and Spectrometer (NIRI) mounted on 
the Gemini-North 8-m telescope in the JHKs bands starting 
on 2007 July 24.566 UT (2.69 hr after the burst); see Ta- 
ble [T] The observations were obtained in excellent seeing 
conditions, « 0.35" in the A'j-band. The data were reduced 
using the gemini package in IRAF, and individual stacks 
were created in each filter Inspection of the images reveals a 
point source coincident with the south-east edge of the puta- 
tive host galaxy; see Figure[T] 

Optical observations were obtained with the Gemini Multi- 
Object Spectrograph (GMOS) mounted on the Gemini-North 
8-m telescope in the gi bands starting on 2007 July 24.551 
UT (2.33 hr after the burst; Table [T]!. The data were reduced 
using the gemini package in IRAF. The combined /-band 
image centered on the location of GRB 070724A is shown in 
Figure |2] An apparent extension is seen in the same location 
as the near-IR source. 

We obtained follow-up near-IR observations of the GRB 
with Persson's Auxiliary Nasmyth Infrared Camera (PANIC) 
mounted on the Magellan/Baade 6.5-m telescope on 2008 
November 18.15 UT in the band (Table[T]i. The individual 
images were dark-subtracted, flat-fielded, and corrected for 
bad pixels and cosmic rays. We then created object masks, 
which were used to construct improved flat fields for a second 
round of reduction. The data were finally registered, shifted, 
and co-added. The resulting combined image is shown in Fig- 
ure [T] and clearly reveals that the point source visible in the 
NIRI images has faded away. 

Similarly, late-time optical gi-hand observations were 
obtained with the Low Dispersion Survey Spectrograph 
(LDSS3) mounted on the Magellan/Clay 6.5-m telescope on 

2008 December 7. 14 UT (Table[T]i. The data were reduced us- 
ing standard procedures in IRAF. The resulting /-band image 
is shown in Figure |2] As in the case of the late-time near-IR 
observations, the extension seen in the early GMOS observa- 



tions is no longer detected. 

To confirm the fading afterglow and to obtain accurate pho- 
tometry and astrometry we perform digital image subtraction 
on the NIRI and PANIC K.-hand images and on the GMOS 
and LDSS3 gi band images with the ISIS package (Alarl 
2000), which accounts for variations in the stellar point- 
spread function (PSF). We adopt the PANIC and LDSS3 im- 
ages as templates with zero afterglow contribution since they 
were obtained about 1 .3 yr after the burst. The resulting resid- 
ual Ks- and /-band images are shown in Figures [T] and |2] re- 
spectively, and clearly demonstrate that the point source co- 
incident with the host galaxy has faded away. We therefore 
conclude that this source is the afterglow of GRB 070724A. 
We additionally performed image subtraction on the two NIRI 
TiTj-band observations and find that the source has faded be- 
tween mid-epochs of 2.832 and 3.696 hr after the burst at a 
4(7 confidence level (Figure[T]l. No residual is detected in the 
subtracted g-band image. 

3.1. Absolute and Differential Astrometry 

We determine the absolute position of the afterglow from 
the NIRI and GMOS residual images using the SExt ractor 
software package^. Since our NIRI images do not contain 
any 2MASS stars, we first perform an astrometric tie of the 
GMOS /-band image relative to USNO-B (using 13 common 
objects with a resulting rms of 0.15") and then tie the NIRI as- 
trometry to the /-band image (using 20 common objects with a 
combined total rms of 0.18"). The optical afterglow is located 
at a = 01^51"14.07P, 5 = -18°35'39.33", while the near-IR 
afterglow position is Q^ = 01^'51™14.066^ 5 = -18°35'39.34" 
(J2000). These positions are consistent within the uncertainty 
of the astrometric tie. The optical/near-IR afterglow is offset 
by about 0.5" and 0.9" relative to the X-ray positions from 
IE vans et all ( |2009|) and [Butler (2007), which have 90% con- 
tainment errors of 1.6" and 1.7", respectively. 

The detection of the optical/near-IR afterglow makes the as- 
sociation of GR B 070724A w ith the previously-proposed host 
galaxy secure (iBergeij|200 9'). and allows us to precisely mea- 
sure the offset between the GRB and host center We perform 
differential astrometry on the NIRI images and find that the 
offset is 0.34 ± 0.01" east and 0.75 ± 0.01" south of the host 
center, corresponding to a radial offset of 0.82 ± 0.01". The 
uncertainty reflects the centroiding accuracy of both the af- 
terglow and host, which we determine using SExtractor. 
At a redshift of z = 0.4571 the scale*" is 5.785 kpc arcsec"', 
and the offset is therefore 4.76 ± 0.06 kpc. This is similar to 
the offset s measured for previous short GRBs with optical af- 
terglows ("Berger et al ."20051: iFox eTani2005l;[Soderberg et alj 
2006; D'Avanzo et al. 2003). The location of the burst in the 
late-time PANIC and LDSS3 images does not exhibit any ex- 
cess emission; see Figure[3] 

3.2. Photometry 

Photometry of the afterglow was performed on all residual 
images using photometric standard stars that were observed 
in conjunction with the PANIC and LDSS3 observations. 
We find that the afterglow had K, = 19.59 ±0.16 mag and 
Ks = 19.64 ±0.17 mag in the first and second NIRI epochs, 
respectively, and / = 23.79 ± 0.07 mag in the GMOS observa- 
tion. The 3(7 limit on the g-band magnitude is g > 23.5 mag, 

^ |http : //sext ractor . source forge ■ net/| 
^ We use the standard cosmological parameters. Ho = 11 km s"' Mpc"', 
n,„ = 0.27, and = 0.73. 
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determined by placing ^10^ random apertures on the resid- 
ual image and using the width of the resulting Gaussian flux 
distribution as Icr. The near-IR magnitudes are quoted in the 
Vega system, while the optical magnitudes are given in the AB 
system. Since the g-band limit is shallower than the detected 
/-band magnitude, it provides no meaningful constraints on 
the properties of the afterglow. 

The observed Ks- and /-band magnitudes correspond to 
fluxes of9.3±1.5 ^lJy, 8.9±1.5 fiJy, and l.liO.l ^iJy, re- 
spectively. We stress that the uncertainty in the flux of the 
near-IR afterglow is dominated by the convolution with the 
PANIC image, which was obtained under worse seeing con- 
ditions than the NIRI images. To assess the statistical uncer- 
tainty in the afterglow flux we note that a stellar point source 
with Ks = 19.62 mag, identical to the afterglow brightness, lo- 
cated near the afterglow position has a Icr uncertainty of 0.03 
mag in the NIRI images. This greater depth in the NIRI data 
allows us to detect a significant fading between the two NIRI 
observations despite an uncertainty of 0. 16 mag relative to the 
PANIC observation. 

We note that the afterglow is also clearly detected in the 
J- and //-band images from NIRI. However, due to the lack 
of late-time template images we cannot robustly measure its 
brightness in these bands. Still, a color-composite image re- 
veals that the afterglow is redder in the near-IR bands than the 
rest of the host galaxy; see Figured 

4. AFTERGLOW/MINI-SUPERNOVA PROPERTIES 

A comparison of our afterglow near-IR flux measurements 
to con temporaneous limits from UKIRT observations by 
iKocevs ki et al. (2009) reveals that the afterglow is about a fac- 
tor of three times fainter than the UKIRT upper limits. These 
authors also find a limit on the optical emission at 0.93 d after 
the burst of /v.; ^ 0.3 /iJy. A comparison to our detected /- 
band flux at 0. 12 d indicates that the afterglow temporal decay 
index is a < -0.6 (/v oc f"), typical of GRB afterglows. 

On the other hand, a comparison of our contemporaneous 
Ks- and /-band fluxes reveals an unusually steep spectral in- 
dex, P = -2.0 ± 0.2 (/v oc v^). Typically we expect /3 ~ -0.6 
to -1 .2 for a wide range of electron power law indices and val- 
ues of the synchrotron cooling frequency ( Sari et al. 1998). 

The unusually red afterglow can be explained in two ways. 
First, the optical emission may be suppressed by extinction 
within the host galaxy. To reconcile the observed fluxes with 
a typical spectral index of /3 « -0.6 requires E{i-Ks)obs ~ 1 
mag, or a rest-frame Ajj?^' « 2 mag for a Milky Way extinction 
curve. Such a large extinction seems unlikely given the loca- 
tion of the afterglow at the edge of the host galaxy. However, 
a large average value of E{B-V) « 1 .2 mag (i.e.. Ay « 4 mag) 
has been inferred for the host galaxy based on its ratio of H7 
and H/3 emission lines (Kocevski et al. 2009), indicating that 
extinction may indeed play a role in suppressing the optical 
emission. 

We further investigate this possibility by comparison to the 
X-ray afterglow brightness at the time of the optical observa- 
tions, /v,x(lkeV) « 0.04 ^Jy (Ziaeepouret al. 2007). This 
leads to an optical to X-ray spectral index of /3o.x ~ -0.5, 
which marginally qual i fies GRB 070724A a s a "dark burst" 
(iJakobssonet al.l l2004t ICenko et al.1 l2009bl) . On flie other 
hand, the near-IR to X-ray spectral index, /3nir,x ~ -0.7, is 
consistent with a typical afterglow. Thus, the comparison of 
the optical/near-IR and X-ray afterglow emission is consistent 
with a standard af terglow origin and sig nificant dust extinc- 
tion. We note that iKocevski et all (l2009h find excess absorp- 



tion in the early X-ray data, but attribute this result to rapid 
variations in the X-ray flux and spectral hardness. In light of 
the possible significant dust extinction, the excess photoelec- 
tric absorption may indeed be real. 

An alternative explanation is that the near-IR flux is dom- 
inated by a different source of emission than the after- 
glow. In particular, in the context of a compact object 
merger, the emission may be due to the decay of radioac- 
tive material synthesized in a sub-relati vistic outflow, the so- 
called Li-Paczynski mini- s upern ova dLi & P aczvriski 19981 
Rosswog & Ramir ez-RuizI (20021 see also Kocevski et aQ 
20091) . In the formulation of Li & P aczynskil fl998). the emis- 
sion from such a mini-SN is described by a peak luminosity 
(Lp): 

Lp^2x lO^V-sMlf (3/3)'/2(k/^j-i/2 g^g s"' , (1) 
a peak effective temperature (Teft p): 

Teff,^ « 2.5 X l0^ffM:l'\3(3T'/\K/Ke)-'^^ K, (2) 
and a peak time (tp): 

fp«lMlf(3/3)-'/2(K/K,)'/2 d (3) 

where / is the fraction of rest mass energy released by the 
radioactivity, M is the ejecta mass in units of M0, /3 = v/c 
is the ejecta velocity, k is the average opacity, ~ 0.2 cm^ 
g"' is the electron scattering opacity, and we use the notation 
X = 10'%. 

For our detected source we use the near-IR luminosity and 
observed time as proxies for Lp and tp, respectively, leading 
to Lp « 10^^* erg s"' and tp « 0.1 d. Using the constraint that 
3/3 < 1 and assuming that k = Ke, we find from Equation|3]that 
M-2 < 10"^ (i.e., M < 10""^ Mq). In conjunction with Equa- 
tion [T] this provides a lower limit of /_3 > 5. The resulting 
lower limit on the effective temperature is Teff,,, > 7 x 10'* K, 
corresponding to a peak in the UV rather than in the near-IR. 
The apparent discrepancy in the spectral peak may be viewed 
as an indication that the observed emission is not due to a 
mini-SN. However, we note that Equations [T][3]correspond to 
the case of a power law decay model with an assumed contri- 
bution from elements with a wide range of decay timescales. 
An exponential decay model, in which a single element dom- 
inates the release of ener gy, may lead to distin ctly different 
luminosity and evolution ( Li & Paczynskil[T998h . 

To summarize, the unusually red optical/near-IR counter- 
part of GRB 070724A, can be explained as a typical afterglow 
with significant dust extinction, Ay°'*' « 2 mag. The alterna- 
tive explanation of a mini-SN leads to an expected peak in the 
UV, but this may suggest that the mini-SN models should be 
revised. 

5. DISCUSSION AND CONCLUSIONS 

Optical afterglow emission has now been detected from 16 
short GRBs, including GRB 070724A. In Figure|5]we plot the 
flux of each optical afterglow at the time of its discovery. For 
GRB 070724A we show both the optical and near-IR fluxes, 
as well as the expected /-band flux extrapolated from the K,- 
band using a typical spectral index of /3 = -0.6. While the 
observed /-band flux is one of the faintest to date, the near- 
IR flux indicates that the afterglow of GRB 070724A is actu- 
ally one of the brightest at the time of its discovery. Indeed, 
only the optical afterglows of GRBs 050724, 060313, 070714, 
and 090510 were brighter, and of these only the afterglow 
of GRB 050724 was discovered on a comparable timescale; 
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the optical afterglows of GRBs 060313, 070714, and 090510 
were all discovered < 20 min after the burst. On a timescale 
of 1 hr to 1 d after the burst, the optical afterglows of short 
GRBs generally have fluxes of ^ 1-10 /iJy, about two or- 
ders of magnitude lo wer than the typical brightness of lon g 
GRB afterglows (e.g, iKann et alJlIOOSt ICenko et al.ir2009bl) . 
From the existing distribution we conclude that the detec- 
tion of a substantial fraction of short GRB afterglows requires 
optical/near-IR observations to otab ^ 24 mag within ~ 0.5 d. 

The unusually red afterglow of GRB 070724A can be ex- 
plained with a substantial rest-frame dust extinction, Ay°'*' « 2 
mag. This value is larger than the typical extinction in- 
ferred for rnost lon g-duration GRBs, Ay ^0.1-1 mag (e.g., 
iPerlev et"ani2009ab . and indeed the optical to X-ray spectral 
index, /3ox ~ -0.5, m arginally qualifies GRB 070724A as a 
dark burst (iJakobsson et al. 2004; Cenko et al. 2009b). Since 
the GRB is located on the edge of its host galaxy, it is likely 
that the extinction arises in the local environment of the burst. 
This implies that the progenitor system was not ejected from 
the host galaxy into the halo or intergalactic medium. Instead, 
the large extinction may point to an explosion site within a star 
forming region, or alternatively that the progenitor system it- 
self produced the dust (for example, a binary system with an 
evolved AGB star). The possibility that some short GRBs are 
obscured by dust has important ramifications for the nature 
of the progenitors, and can also serve to localize the bursts to 
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TABLE 1 

Log of Near-IR and Optical Observations of GRB 070724A 



Date 


A? 


IClCaCUpc 


Tti c tni tn f»Ti t 


Filter 


P YT^r\ci It*** c 


^FWHM 


mag 


f 


(UT) 


(d) 








(s) 


(") 




(MJy) 


2007 July 24.549 


0.094 


Gemini-N 


GMOS 


8 


2x 180 


0.71 


>23.5 


< 1.5 


2007 July 24.551 


0.097 


Gemini-N 


GMOS 


i 


2x 180 


0.53 


23 .79 ±0.07 


l.litO.l 


2007 July 24.572 


0.118 


Gemini-N 


NIRI 


Ks 


15 X 60 


0.35 


19.59±0.I6 


9.3 ±1.5 


2007 July 24.585 


0.131 


Gemini-N 


NIRI 


J 


15 X 60 


0.45 


a 


... ^ 


2007 July 24.596 


0.142 


Gemini-N 


NIRI 


H 


15 X 30 


0.46 


a 


... ^ 


2007 July 24.608 


0.154 


Gemini-N 


NIRI 


Ks 


15 X 60 


0.35 


19.64±0.17 


8.9 ±1.5 


2008 Nov 18.17 


481.6 


Magellan 


PANIC 


Ks 


54 X 20 


0.47 


* 


... b 


2008 Dec 7.13 


500.6 


Magellan 


LDSS3 


g 


2x240 


0.94 


b 


... b 


2008 Dec 7.14 


500.6 


Magellan 


LDSS3 


i 


3 X 120 


0.58 


... b 


b 



Note. — "No templates are available in the J and H bands and as a result we cannot measiu-e the flux of the 
afterglow. 

* The flux of the afterglow in the R\NIC and LDSS3 images is assumed to be zero. 
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Fig. 1 . — Near-IR i^-band images of GRB 070724A. In each row we display the afterglow image, corresponding template image, and residual image from 
ISIS. The near-IR afterglow is clearly visible in the residual images relative to the final PANIC epoch, as well as in the subtraction of the two NIRI epochs. 
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Fig. 2.— Optical i-band images of GRB 070724A. 
visible in the residual image. 
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Fig. 3. — Near-IR A:,-band images from NIRI and PANIC of the host and afterglow of GRB 070724A. The cross-hairs on the PANIC image mark the location 
of the afterglow and indicates that GRB 070724A occurred on the outskirts of the host galaxy. The offset measured from the NIRI image is 4.76 ± 0.06 kpc. 
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Gemini/NIRI 
J/H/Ks 
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Fig. 4. — Color composite image of tiie field of GRB 070724A using / (blue), H (green) and Ks (red) images from NIRI. Although we cannot clearly measure 
the afterglow flux in the J and H bands due to the lack of late-time template images, we find that the afterglow is clearly redder than the host galaxy. 
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Fig. 5. — Optical flux at the time of discovery for all 16 short GRBs with optical and near-IR afterglows (including GRB070724A). The solid stars mark the 
I- and ^Tj-band fluxes of the afterglow of GRB 070724A, while the open star is the optical flux in the /-band extrapolated from the near-IR with a spectrum of 
iv oc v^'^ , typical of GRB afterglows. Data for other short GRBs are taken from the literature (Hjorth et al. 20()5b 
Se Uaarte Postiso et al. 2006; Levan et al. 2006; Romins et al. 2006; D'Avanzo et al. 2009; Stratta et al. 2007 
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